Coulomb Interaction and the Fermi Liquid State: Solution by Bosonization 
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We investigate the effects of the Coulomb two-body interaction on Fermi liquids via bosonization. 
The Coulomb interaction is singular in the limit of low momentum transfer, and recent interest in 
the possibility that some singular interactions might destroy the Fermi liquid state motivate us to 
reexamine it. We calculate the exact boson correlation function to show that the Fermi liquid state is 
retained in the case of Coulomb interactions. Spin and charge degrees of freedom propagate together 
at the same velocity and collective charge excitations (plasmons) exhibit the expected energy gap 
in three dimensions. Non-Fermi liquid behavior occurs, however, for a super long-range interaction 
studied recently by Bares and Wen. 
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I. INTRODUCTION 



> 

. The many body problem of fermions interacting via the Coulomb interaction is central to condensed matter physics. 
The standard method for studying this problem is the summation of ring diagrams which arise in perturbation theory: 
the random phase approximation_(RPA)ld. Recently, another approach to this problem based on the renormalization 
group (RG) has been developed 2 ]. Both the RPA and the RG calculations assume that the Fermi quasiparticle 
propagator retains the Fermi liquid form with a simple pole. This assumption is shown to be self-consistent since 
in both approaches the bare Coulomb interaction is screened down to a short-range form. On the other hand, it is 
desirable to develop tools which do not rely on the intermediate assumption of a Fermi liquid propagator to ascertain 
whether non- Fermi liquid states are solutions of.sy.stems with singular interactions. 

Bosonization in dimensions greater than onetTu is particularly well suited to the study of singular long-range 
. interactions as the realm in which it is applicable, for low-energy excitations near the Fermi surface, is precisely 
where these interactions are strongest. Furthermore bosonization does not rely upon a Fermi liquid form for the 
quasiparticle propagator; for example, it encompasses the possibility of spin-charge separation. Our results, however, 
O ■ are unambiguous: we find that the Fermi liquid state is the solution to the Coulomb problem. By using the bosonization 
transformation to determine the fermion quasiparticle propagator, we obtain well-known results for the fermion self 
energy: the imaginary part is proportional to w 2 ln|tj| in two dimensions and just uj 2 in three dimensions. We 
emphasize that the bosonization method yields non-perturbative information, so a natural next step would be to use 
^ . it to study the effects of transverse gauge interactions. 



II. BOSONIZATION OF THE TWO-BODY INTERACTION 



We begin with the bare Hamiltonian for electrons or other fermions interacting via the two-body Coulomb or other 
long-range longitudinal interaction V(q): 

H = zZ Ck ^ C k« + y ( q ) ^ C P+q C P/3 Ck +q« co 

k k,p,q 

where = k 2 /2m for the case of a Fermi gas and there is an implicit sum over repeated spin indices.^ and (3. Then 
we integrate out the high-energy Fermi degrees of freedom with the use of the renormalization groupBH The resulting 
low-energy effective theory is expressed in terms of quasiparticles ipk a which obey canonical anticommutation relations 
and which are related to the bare electron operators via the wavefunction renormalization factor Z^: 

V'ka = Z k 1/2 c ka (2) 
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for momenta k restricted to a narrow shell of thickness A around the Fermi surface: kp — A/2 < |k| < kp + A/2. The 
resulting low-energy Fermi liquid Hamiltonian is: 

h fl = e ~k ^ Q + ^ J ( s ; -<0 J ( T ; <*) + { re § ular terms > ( 3 ) 

k ° S,T;q 

Here ?k = k 2 /2m* incorporates the mass renormalization. The charge current in a given Fermi surface patch S [where 
S = (9 1 cj>) in three dimensions] is defined by: 

J(S; q) = 9 ( S ; k + q) 9 ( S ' k ) {^k+ q 1><* - «k} • (4) 

k 

Here 0(S;k) — 1 if k lies inside a squat box of dimensions A x A D_1 centered at S and equals zero otherwise. The 
regular terms in Eq. (||) do not diverge in the q — > limit and consist of both large-g exchange processes and effective 
interactions generated by the high-energy electrons that have been integrated out. We assume in the following that 
no superconducting or charge or spin density wave instabilities arise from the regular terms. nn 

We now review the main aspects of bosonization and refer the reader to two previous papersO'B for more details. In 
D-spatial dimensions, the Fermi quasiparticle fields of spin-er, f/v, may be expressedta in terms of the abelian boson 
fields 6 a as: 



</v(S; x) = -±= J -e lk - x exp{^^(S; x)} O(S), (5) 

where the dependence on time, t, is included implicitly in the spatial coordinates x and S labels the patch on the 
Fermi surface with momentum kg. Vol is the volume of the system which equals L D in D-dimensions; the factor of 
y/Vol is introduced to keep the fermion anticommutation relations canonical. Both the ip and (f> fields live inside the 
squat box centered on S with height A in the radial (energy) direction and area A 0-1 along the Fermi surface. These 
two scales must be small in the following sense: kp >> A >> A. We satisfy these limits by setting A = kp/N and 
A = kp/N a where < a < 1 and N — ► oo. The quantity a in the bosonization formula Eq. (|J) is a real-space cutoff 
given by a = 1/A. Here i7 = A d ~ 1 LLJ2tt) d equals the number of states in the squat box divided by A. Finally, O(S) 
is an ordering operator introducedau to maintain Fermi statistics in the angular direction along the Fermi surface. 
(Anticommuting statistics are obeyed automatically in the radial direction.) 

With the connection Eq. (|J) between the fermion and boson fields we may bosonize the free Hamiltonian. The 
result is quadratic in the <j> fields: 

H = vp Y J d D x^ a (S;x) j^Z-fc F j ^ Q (S;x) 

J2 fd D x{(h s -W) Q (S;x)} 2 . (6) 



s 

2-7T Vp 

n Vol 



The total bosonized Hamiltonian may be written as H = H c + H s , .exhibiting the factorization into charge and spin 
sectors. The charge Hamiltonian is bilinear in the current operatorsu J(S; q): 

H - = \ E E ^( s > T ; <*) J ( s ; J ( T ; <0 • ( ? ) 

S,T q 

Long-range interactions V(q) are incorporated into V c as matrix elements that couple currents in different patches: 

V C (S, T; q) = ± fi" V S^ 1 + ^ V(q) . (8) 



The charge currents are related to the-charge boson field <f> = <f>^ + <f>i by J(S; x) = v47rns ■ V0(S; x). They obey the 
equal-time U(l) Kac-Moody relations^: 

[J(S; q) , J(T; p)] = 2 Q S^ 1 ^ +p>0 q • n s . (9) 

Similarly, the spin-sector is described by a Hamiltonian which has a free part: 
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H °=2 EE^^^-q) (io) 

S q 

and an interaction part that depends on the regular Fermi liquid spin-spin coefficients f s . Here the abelian spin 
currents J z commute with the charge currents and are expressed in terms of the spin boson field (f> z = 4>-\ — 4>l by 
J Z (S; x) = \/47rns • V0 z (S; x). The interaction term has two parts: a term that couples the z-component of the spin 
currents in different patches and a term that couples the x- and y- components which has the form: 

Y, /.(S,T) f d D x cos{^p [0,(S;x)-0,(T;x)]} . (11) 

S,T J 

We note that the coefficients f s are invariant under the renormalizatipn group transformations despite the fact that 
the above term resembles the form of the bosonized BCS interaction^. The crucial difference is that here the boson 
fields do not appear in pairs at opposite points of the Fermi surface. The /3-function for the coefficients f s therefore 
equals zero. For simplicity we set / S (S, T) = in the following. 



III. QUANTIZED BOSONS 

In this section we calculate the exact boson Green's function in the charge sector. To simplify the calculation we 
restrict our attention to the case of spherical (circular in two dimensions) Fermi surfaces and set the Fermi velocity 
equal to one. None of these simplifications is essential. First we write the charge currents in terms of boson operators 
that satisfy canonical commutation relations. The choice: 



J(S; q) = y/2Sl |n s -q| [o(S; q) 9(n S • q) + a f (S; -q) 0(-n s • q)] (12) 

with 

[a(S;q), ot(T;p)]=^ 1 , (13) 

and 9(x) = 1 if x > and is zero otherwise, satisfies the U(l) Kac-Moody commutation relation Eq. (^). For 
convenience we denote a(S;q) and a(S; — q) by afl(S;q) and ot(S;q), respectively the right and left moving fields. 
The momentum-frequency space propagator 

%Gi(S; q, w) = (ai(S; q, w) aj(S; q, w)) (14) 
is related to the propagator of the <f> fields by: 

(^(S;q,«) ^(S;-cb-w)) = — " (ai(S;q,o;) oJ(S;q,w)) . (15) 

47rns • q 

As the interaction Wq) has no dependence on the patch indices S and T, a straightforward generalization of our 
previous calculationu for the short-range interaction Fq allows us to construct the boson propagator. The exact 
solution of the Dyson equation for the irreducible self energy is given by: 



2A - 1 n s • q V{q) 
The Lindhard function x° m two dimensions is given by: 



*(S;^)= (27r) . i + vw (16) 



q 2n cos(cf>) — x — if] sgn(cj) 
, s f Mx 1 - 1) , ,6>(1 -x 2 ) I 

where x = ^ and N(0) is the density of states at the Fermi energy which in two dimensions equals kp/ir in units 
where vf = 1- In three dimensions the Lindhard function equals: 
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X°(aO=iV(0) |l-|ln(£±l)| ; D = 3. (18) 

Despite the fact that a perturbative expansion!! has been used as an intermediate step to obtain Eq. (0), all terms 
in the expansion have been summed yielding the exact non-perturbative result valid for arbitrary dimension D. For 
example, when the interaction is short-ranged, in D = 1 the resummcd expansion yields the well known exact result 
for a Luttinger liquidcl. 

In the N — > oo limit in which q — ► the self energy simplifies because V^(q) — > oo. For finite- x we obtain: 

S*(S;q,c,) = ^g^ . (19) 

Screening is apparent at finite-x, even in the boson propagator, as the self energy in the small-q limit no longer 
depends on V(q). The velocity is renormalized slightly from its bare value of unity, v' F = 1 + 0(-^-) D ^ 1 , and the 
boson lifetime is now finite because of scattering into different patches. Note, however, that these changes represent 
irrelevant corrections as the self energy Eq. ([16]) scales to zero as A — > 0. In particular, the pole in the boson 
propagator remains unchanged as N — > oo. In the opposite w-limitQ of x — > oo the denominator in Eq. ( |l6| ) vanishes 
at frequencies corresponding to the eigenvalues of the collective mode equation; consequently the self energy diverges 
at these frequencies and in the three dimensional Coulomb case the plasmons acquire a gap as expected!! 

IV. FERMION QUASIPARTICLE PROPERTIES 



In the previous section we saw that long-range longitudinal interactions modify the boson propagator. Though 
BCS scattering processes were ignored, small-angle scattering processes made the boson lifetime finite. With these 
results we can use the bosonization formula Eq. (j^) to infer the fermion quasiparticle self energy. Since bosonization 
is carried out in (x, t) space we must carry out three operations. First we Fourier transform the charge and spin boson 
propagators into real space and use the abelian relation </>j,j = |(<^>± <fi z ) to construct the boson propagator for, say, 
up spins. Next the exponential of the resulting expression yields the fermion propagator in real space. Finally an 
inverse Fourier transform of the fermion propagator back into momentum space allows us to extract the self energy. 

It is difficult technically to perform these steps in all generality. It will be sufficient for our purposes to estimate 
the imaginary part of the fermion self energy. The charge boson Green's function may be written in configuration 
space as: 



G (S;x,t) ee (0(S;x,t) 0(S;O,O) - 2 (S; 0, 0)) = i FT {G^} 

^(o) 

V 1 ^<t> 



l FT{Gf+G^} (20) 



where FT represents the Fourier transform operation that converts the variables (q, u>) to (x, t). In the second line, 
iGfy which is given by Eq. (|l5|), has been expanded in powers of the vanishingly small parameter A/kp. The Fourier 



transform of the real term, i FT{g£ }, yields the free propagator: 

O 2 in 

(<XT; x) 0(T; 0) - 2 (T; 0)) = — ln( ,° ) ; |x ± A| « 1 

47r x ■ riT — t 

-> -oo ; |x_lA| » 1 (21) 

where x^ denotes spatial directions perpendicular to the surface normal fix. The contribution to the imaginary part 
of the fermion self energy comes from G^ 1 ' in the limit of \x\ < 1 for which quasiparticle damping occurs: 

G^(S;q,w) = - [to - n s ■ q + ir] sgn(^)] 2 

2A D ~ 1 

= ~ A /r, yd r [w-fis-q + ^sgn^)]- 2 (22) 

47T [2-k) u x (q, w) 

where the second line holds only for finite- a: and small q. The imaginary part of the fermion self energy at the pole 
may now be estimated by-, exponentiating Eq. ( po|) and expanding to first order in the imaginary part. Following the 
steps in a previous papeia we obtain for D = 2: 
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which form we recognize immediately from previous work on the quasiparticle lifetime within the RPA approximation^. 
This quantity is always negative since \lo\ << A. In three dimensions the calculation yields: 

Im S^MUic 4 w2 + 0(u 3 /k 2 F ); (24) 
Hp 

the appearance of the momentum cutofLA in this expression coincides with that found in traditional Fermi liquid 
theoryE3. In the classic RPA calculation^ the cutoff A is replaced by an energy of order the plasma frequency due to 
a different weighting of the high-energy states. 

In the case of Coulomb interactions the weight of the quasiparticle pole, Zp , remains non-zero at the Ferrpi, surface. 
In contrast, super long-range interactions can destroy the quasiparticle pole. For example, Bares and WentHl studied 
an interaction in two dimensions with a logarithmic potential given in momentum space by V(q) = g/q 2 and showed 
that Zp = within RPA. Note that g is a momentum scale of order the Fermi momentum and the plasmon gap is 
non-zero and equals y / gkp/2TT due to the super long-range nature of the interaction. In fact, this system may form 
a Wigner crystal; we ignore any such 2kp instabilities in the following analysis. We confirm that the quasiparticle 
pole is destroyed in the super long-range case via bosonization. One might be tempted to compute Zp by using the 
Kramers-Kronig relation to derive the real part of the self energy from the imaginary part estimated above. This 
procedure however, is unreliable, as the Kramers-Kronig relations involve an integral over the entire energy range 
whereas the calculation breaks down at energies greater than A. Instead we directly calculate the real-space Fermi 
two-point function. To do this we compute the Fourier transform of Eq. ( p2| ) in the |i| > 1 region for which the self 
energy Eq. ( |l6| ) is purely real. For simplicity we also expand the Lindhard function in powers of 1/x: 

xV) = -^r + o(i/* 4 ). (25) 

The important point is that, in the case of the equal-time two-point function, the integral over q diverges logarithmi- 
cally in the infrared due to a factor of q 2 appearing in the denominator of the momentum integral. Setting t = we 
obtain: 

, si /o ■> , /- \ f d 2 q exp(iq-x) f x duo 1 1 
lnG/(S;x) « -ln(n s • x) + / — ^ '- / — — — — — —-^ (26) 



9c <| q |<A ( 27r ) 2 I 2 J\q\ * Vff- k F /(2Truj 2 ) [u - n s • q + ir\ sgn(w)] 2 

Since to > |q| in the integral, we make the approximation of replacing [u> — ns ■ q + if] sgn(w)] 2 with oj 2 + ir( in the 
denominator and carry out the integrations. Upon exponentiating the boson Green's function we obtain the Fermi 
quasiparticle propagator in configuration space: 

G/(S;x) a e^' 1 ; |n s • x|A » 1, |x_l|A « 1 (27) 



n s 



where 

The appearance of the infrared cutoff q c = 0(1/ L) in the numerator reflects the super long-range nature of the 
Bares and Wen interaction. In a Luttinger liquid with short-range interactions, a length scale of order the lattice 
cutoff would appear in place of the system size L. The appearance of the anomalous exponent £ is a consequence 
of bosonization which treats the interaction japn-perturbatively and thus improves upon RPA which gives instead 
a logarithmic dependence on the system sizaia. Evidently, the quasiparticle pole has been destroyed by the super 
long-range interaction. An evaluation of the fermion propagator in the opposite limit of x = and at large times t, 
however, yields the usual i _1 Fermi liquid form as screening is effective at low frequencies. Likewise, thermodynamic 
properties which are defined in this equilibrium g-limit such as the specific heat also have the usual Fermi liquid form. 
Note that in either limit the propagator is odd (apart from the phase factor) under (x, <)—»•— (x, t) as demanded by 
Fermi statistics. 

Repeating the calculation in the case of Coulomb interactions in either two or three dimensions we easily find that 
the propagator has the standard Fermi liquid form. 
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V. CONCLUSION 



We examined the effects of long-range interactions on fermion liquids by bosonization in dimensions greater than one. 
We find that the Fermi liquid state is the only solution to the problem of a degenerate gas of fermions interacting via the 
Coulomb two-body interaction in two and three spatial dimensions. Bosonization allows us to go beyond an assumed 
Fermi liquid form for the quasiparticle propagator. Indeed, the fact that the bosonized Hamiltonian separates into 
charge and spin parts, H = H c + H s , leads to the possibility that, as in one dimensionta, the quasiparticle propagator 
might also exhibit spin-charge separation, especially in the case of the Coulomb interaction which is singular. Spin- 
charge separation in dimensions larger than one would, however, destroy the Fermi liquid as the key element, the 
existence of a pole in the single-particle Green's function with non-zero spectral weight, would be replaced by a branch 
cut and the pole would be destroyed. This does not happen because, as we saw at the end of section III, the location 
of the pole of the boson propagator is unchanged from its free value in the A — > limit. Consequently the spin and 
charge velocities are equal and both degrees of freedom propagate together in the usual quasiparticle form. On the 
other hand, the Fermi liquid form is destroyed in the case of the super long-range interaction studied by Bares and 
Wen. 
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